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ABSTRACT 



Experimental studies were conducted on unfilled visco- 
elastic materials to evaluate the fracture postulate which 
states that crack velocity is a unique function of stress 
intensity factor. To do this, fracture characterization 
tests were conducted for constant load, constant displace- 
ment, and constant displacement rate histories. A cyclic 
study, a more sophisticated application of the constant dis 
placement rate case, was also conducted as a further evalua 
tion of the fracture postulate. 

The analytical theory used in obtaining viscoelastic 
stress intensity factors was based upon linear theory and 
sharp crack geometries. Moderate correlation was obtained 
among all test cases, notwithstanding the presence of 
significant specimen strains and finite crack tip radii. 
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I 



INTRODUCTION 



Over the years, fracture studies have led to the devel- 
opment of relationships among geometry, loading and crack 
length that prevail at the incidence of rupture and the 
onset of crack propagation. It has been found useful to 
define a stress intensity factor, , which increases as 
the load on the specimen increases until a critical value, 

is reached* At this point the crack rapidly propagates, 
and it has been found that is a material property, such 

that in every structure when this condition is reached, the 
crack becomes unstable. 

In viscoelastic materials, the crack is found to propa- 
gate for almost any load level. The velocities are infini- 
tesimally small for low loads and increase with increasing 
load. The question then becomes: at what velocity does 

the crack propagate for any given load, not at what load 
level does rapid crack propagation occur. A stress inten- 

. V 

sity factor for the viscoelastic case, , can be obtained 

by the correspondence principle, where the stress intensity 
factor, , for the elastic case in the Laplace transform 

domain is inverted into the time domain to give the visco- 

V 

elastic stress intensity factor, contains information 

relative to the specimen geometry, material, and boundary 
conditions . 



7 



The fracture postulate for viscoelastic materials can 



be stated in the following way: 



If 



tKi [K^ ^’^^'bODY 2 



then [4 (t)lgo„y 3 _ - (4 (t))ggpy 2 



Hertzler [Ref. 2] conducted constant stress and constant 
displacement tests on filled viscoelastic materials in an 
attempt to substantiate this relation. He developed frac- 
ture characterizations by plotting stress intensity factors 
versus crack velocities for each test. His results showed 
the constant displacement and constant stress fracture 
characterizations diverging with increasing stress intensity 
factor. This disparity was attributed to the dewetting 
behavior of the material. 

This study makes a further attempt to prove the fracture 
postulate using an unfilled viscoelastic material, since 
unfilled materials do not experience dewetting. Fracture 
characterizations were developed, using the same procedure 
as Hertzler's, from three experimental tests: the constant 

load case, the constant displacement case, and the constant 
displacement rate case. A cycling application of the con- 
stant displacement rate test was also performed as an 
added history common to viscoelastic analysis and as a fur- 
ther test of the fracture postulate. These test results 
were then compared, to evaluate the fracture postulate. 
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II . EXPERIMENTAL EQUIPMENT 



A, TEST SPECIMENS 

The unfilled viscoelastic material used in testing was 
a carboxy terminated^ poly butadiene (CTPB) polymer composed 
of 94.97% Butarez and 5.03% Hx868. The samples were fabri- 
cated by the Naval Weapons Center, China Lake, California. 

The specimens were prepared for testing by cutting a 
three-fourths inch horizontal crack in the sample as shown 
in Figure 1. Specimens could then be used for crack lengths 
up to three inches with less than 1% difference between the 
far-field stress in the undisturbed region and the average 
stress based on net cross section [Ref. 1] . 
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Figure 1. Geometry of Test Specimens 
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B. LABORATORY HARDWARE 



Experimental testing was done with a table model 
Universal Instron. A temperature chamber was built to 
fit the Instron so that specimen temperature could be 
regulated during the test (see Figure 2). The Instron 
was equipped to directly conduct constant displacement, 
constant displacement rate, and cyclic tests. 




Figure 2. Universal Instron with Temperature Chamber. 



For the constant load test, the apparatus was adapted 
(Figure 3) so a constant load could be suspended beneath 
the temperature chamber. Since this case required a 
uniform stress state on the loaded boundary, the apparatus 
was built so that the center of gravity of the applied 
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load could be maintained under the center of the uncracked 



length of the specimen. This provided a constant strain 
along the length of the specimen, insuring a uniform stress 
boundary condition. 




Figure 3. Universal Instron and Temperature Chamber 
Adapted for Constant Load Test 



Specimen temperature was monitored with five thermis- 
tors each located one-eighth inch behind the sample and 
spaced so as to provide a thorough coverage of the 
material surface. Their locations are shown in the photo- 
graph of Figure 4. 

Even though the temperature coils were exposed and 
located relatively close to the specimen ends, minimal 
temperature gradients occurred in the area of the test 
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Figure 4. Temperature Chamber. 

sample. The uniform temperature in the horizontal direc- 
tion was attributed to the large metal fixture in the 
center of the chamber which acted as a heat sink, provided 
much thermal inertia, and tended to stabilize the tempera- 
ture in its immediate vicinityo The temperature chamber 
had a substantial vertical temperature gradient, but this 
feature was not critical since the specimen dimension in 
this direction was small. 

Temperature chamber statistics are given in Table 1. 
The "high*' and "low" data represented the test specimen 
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environment at the maximum and minimum temperatures with- 



in the chamber during one cycle. This high- low-high 
chamber cycle took about five minutes for the 15, 30 and 

50 degree temperatures, while the compressor ran for a 
longer portion of the cycle as the temperature decreased. 
At 70 degrees, the chamber was basically maintaining room 
temperature but providing a chamber environment. At zero 
degrees the compressor was continually running providing 
a relatively constant temperature. 
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Ill 



FJXPERIMENTAL RESULTS 



A. CONSTANT LOAD CASE 

The viscoelastic stress intensity factor (opening 
mode) for the constant load case was found by Lindsey 
[ Ref . 1 ] to be : 



K 



V 

1 



(t) 



V 



3b(t) 

4TT 



F 

A (t) 



( 1 ) 



This test consisted of placing a step load on the 
sample^ which was then maintained constant throughout 
the test. Load levels required to produce crack propa- 
gation created specimen strain levels from 20% to 40%. 

To provide a uniform stress on the loaded specimen 
boundary and to avoid introducing a moment at the crack 
tip as the crack projDagated^ the center of gravity of 
the applied load was continually moved so as alv/ays to 
be under the center of the uncracked length of the sped 
men. Slow crack velocities made this method tractable. 

Since the specimen was finite in length, an increase 
in crack length resulted in a decrease in A(t) and an 
associated increase in material strain. This increased 
the sample width and, by Poisson's effect, reduced the 
sample thickness. V7ith increasing crack length, then, 
the stress intensity factor increased, resulting in a 
history of stress intensity factors for one applied load 
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For the constant load case, the increase in crack 



velocity with an increase in crack length was small. 
Therefore, the crack velocity was experimentally observed, 
and averaged about four times, during each test. The 
stress intensity factor was computed based on the average 
crack length for a particular crack velocity. Stress in- 
tensity factors versus crack velocities were then plotted 
for specific testing temperatures. Since the visco- 
elastic material was thermor eological ly simple [Ref. 2] , 
each temperature curve was ultimately shifted into a 
smooth master fracture characterization curve. 

Individual temperature test results are given in 
Figure 5, and the master fracture characterization curve 
is given in Figure 6. Test data had relatively little 
scatter, and the experimental temperature shift factors 
had good correlation with relaxation modulus shift 
factors [see Figure 6] . 

B. CONSTANT DISPLACEMENT CASE 

The viscoelastic stress intensity factor (opening 
mode) for a constant displacement case was found by 
Lindsey [Ref. 1] to be: 



To obtain a constant displacement boundary condition, 
the specimen boundary was displaced at a rate of two 
inches per minute until reaching the desired displacement. 
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Figure 5. Fracture Characterization Curves for the 
Constant Load Case. 
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Figure 6. Fracture Characterization Master Curve for the: Constant Load Case 



I 



V^. Then, by waiting ten times the loading time, the ramp 
tests simulated step displacement loading for which the 
stress intensity factor reduced to Equation (2). As time 
increased after loading, the stress intensity factor and 
corresponding crack velocity decreased as a result of 
relaxation. This yielded a history of fracture* data from 
one test. 

As in the constant load case, the crack velocity was 
experimentally observed and averaged over each time incre- 
ment. Stress intensity factors were computed based on 
the average relaxation modulus for a particular crack 
velocity. Stress intensity factors versus crack velocities 
were plotted as a function of temperature, and finally 
shifted to a master fracture characterization curve. 

Individual temperature test results are given in 
Figure 7 and Figure 8, and the master fracture characteri- 
zation curve is given in Figure 9. Test data had a 
moderate degree of scatter, and experimental temperature 
shift factors had some disagreement with relaxation 
modulus shift factors [see Figure 9] . 

C. CONSTANT DISPLACEMENT RATE CASE 

The viscoelastic stress intensity factor (opening 
mode) for a constant displacement rate case was found by 
Lindsey [Ref. 1] to be: 
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Fracture Characterization Curves for the 
Constant Displacement Case 
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Figure 9. Fracture Characterization Master Curve for the Constant Displacement Case 



where were determined by the procedure given in 

ICRPG Solid Propellant Mechanical Behavior Manual [Ref. 3] 

The constant displacement rate test utilized the 
Instron in its most straight-forward application. One dis 
placement rate was selected for the Instron crossarm, and 
this rate was maintained throughout the complete test. 

This test differed from the constant load test and the con 
stant displacement test in two ways. First, one test 
generated a wide range of crack velocities, prohibiting 
the averaging of crack velocity during the test. Secondly 
this extensive history of data from one test resulted in 
only three different constant displacement rates at one 
temperature being required for sufficient fracture charac- 
terization . 

To accurately represent experimental data, crack 
lengths versus time were plotted. Fracture characteriza- 
tion was then done in two ways. First, crack velocities 
were graphically obtained by measuring the slope of the 
crack length versus time plots, and corresponding stress 
intensity factors were computed. This method permitted 
a direct comparison with previous results. 

Secondly, recognizing the inaccuracy inherent in 
graphically obtaining crack velocities, an incremental 
technique was employed. The crack le ng th- ver su s - time plot 
was divided into discrete time increments. Beginning at 
the first time increment, a stress intensity factor was 
computed from Equation (3), and a corresponding crack 
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velocity was obtained from the master fracture characteri- 
zation curves. The current crack length was then obtained 
from the crack velocity by assuming it to be constant over 
the time increment. This process was continued and the 
resulting crack length-time history was compared with the 
experimental crack length-time history. 

Experimental results based on graphically obtaining 
the crack velocity are given in Figure 10, labeled ''Con- 
stant Displacement Rate." Crack length versus time plots 
are given for the three constant displacement rate tests 
in Figure 11, Figure 12, and Figure 13, and are labeled 
"Experimental, Constant Displacement Rate." A smooth, 
continuous crack lengt h-ver su s- t ime curve resulted from 
the constant displacement rate testo 



D. CYCLIC CASE 



The basic equation for a viscoelastic stress intensity 
factor (opening mode) for the cyclic case v;as given by 
Lindsey [Ref. 1] as: 



(t) = 
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- a . t-T ) 
e 1 






(4b) 



The cyclic case was an extension of the constant dis- 
placement rate case as one crossarm rate was maintained 
throughout the complete test. The direction of crossarm 
travel was reversed, however, when a selected displacement 
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Summary of Fracture Characterization Master Curves 
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Figure 11. Crack Length versus Time for the Constant Displacement Rate Case 
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Figure 12. Crack Length versus Time for the Constant Displacement Rate Case 
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was reached. This resulted in a linear, piecewise con- 
tinuous, even function which could be approximated by a 
Fourier series. 




Figure 14. Basic Cycle for the Cyclic Test. 
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Substituting Equation (5b) and Equation (4b) into 
Equation (4a) resulted in: 



K ''(t) = — 

1 ^ ' TTL t 3b (t) -rr 
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Equation (6) was solved on the computer using 14 terms in 
the Fourier series. 

Experimental data were plotted as crack length versus 
time. This was then compared with the constant load case 
and constant displacement case using the incremental pro- 
cedure described in the constant displacement rate case. 

A crack leng th-ver sus- t ime plot is given for a cycling 
test in Figure 15 (also shown in Figure 16) , labeled, 
”Exper imental , Cyclic Case.” This test produced a stair- 
step crack length history. 
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Figure 15. Crack Length versus Time for the Cyclic Case. 
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Figure 16. Crack Length versus Time for the Cyclic Case 
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IV. DISCUSSION OF RESULTS 



The constant load, constant displacement, constant dis- 
placement rate, and cyclic tests were diverse in nature, 
each having several positive aspects and some having minor 
deficiencies. Each of these four tests probed different 
aspects of material behavior, and collectively they provided 
a thorough and comprehensive examination of the fracture 
postulate . 

The constant load case had two positive features. The 
test occurred over a relatively long period of time, simpli- 
fying laboratory techniques and stabilizing test conditions, 
and the theoretical stress intensity factors were indepen- 
dent of external relaxation data. Detracting from this 
test, however, was the mechanical apparatus used to apply 
a constant load. Unless carefully manipulated, it would 
induce a moment at the crack tip. 

The constant displacement case had several attractive 
features. Stress intensity factors were calculated using 
relaxation data obtained from constant displacement tests. 
Since the test material was an elastomer, the reduction in 
crack velocity due to relaxation occurred over a substan- 
tially longer time than a similar test using solid propel- 
lant. Hertzler [Ref. 2] found this test virtually impos- 
sible to use for propellants because rapid relaxation 
caused the crack velocities to go to zero before data 



could be taken. Finally, the Instron was capable of a 
consistent application of the constant displacement case. 

The constant displacement rate case utilized the 
Instron in its most fundamental mode of operation. Since 
few tests were required for fracture characterization, 
these tests were repeated. Excellent consistency v/as 
obtained between like tests. A primary disadvantage was 
that the stress intensity factor calculation required 
relaxation data. The only relaxation data available were 
those obtained from a constant strain case. 

The comments pertaining to the constant displacement 
rate case were relevant to the cyclic case. The cyclic 
case was a more thorough test, however, as disparities 
between experimental and theoretical results were accumu- 
lated over each cycle, and a discrepancy not apparent in 
one cycle (equivalent to one constant displacement rate 
test) was apparent after several cycles. 

The summary of the master fracture characterization 
curves given in Figure 10 graphically compares the con- 
stant displacement, constant load, and constant displace- 
ment rate cases. This shows a correlation of results 
among the three tests within an error band of % measured 
on the vertical. Figure 11, Figure 12, and Figure 13 com- 
pares the three constant displacement rate cases with the 
constant load case and the constant displacement case 
using the incremental technique. These results were con- 
sistent with Figure 10, as moderate correlation was 



33 



obtained in all test cases. Figure 15 and Figure 16 com- 
pare the cyclic case with the constant load case and the 
constant displacement case, respectively, using the incre- 
mental technique. The cyclic test had good correlation 
with the constant load test through four cycles. It would 
be difficult to maintain this agreement through a large 
number of cycles, however, as the incremental process 
resulted in the error in one cycle being added to succeed- 
ing cycles. This manifestation of the cyclic test is 
illustrated in Figure 16 where the constant displacement 
and constant displacement rate cases became further sepa- 
rated with increasing cycles, because only a moderate 
agreement existed between the two test cases in any one 
cycle. 

The analytical theory used in obtaining viscoelastic 
stress intensity factors was based upon linear theory in 
that the change in the sample dimensions, when it was 
loaded, was assumed to be small. In actual test cases, 
however, specimens experienced from 20% to 40% strain, 
causing considerable variation in s amp le width and thick- 
ness. This represented a significant departure from the 
linear case. 

A second departure from the analytical theory occurred 
when the crack developed a finite crack tip radius. The 
classical theory for the stress intensity factor was 
developed for an infinitely sharp crack. When the unfilled 
viscoelastic test specimens were loaded, however, the crack 
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tip was approximately elliptic* Figure 17 shows two 
specimens at different strains under typical loading 
situation. These photographs indicate not only a finite 
crack tip radius but also a change in radius with a 
corresponding change in material strain. Since the con- 
cept of a stress intensity factor as a function of crack 
tip radius has not been developed, one can only specu- 
late on its effect. 



35 



jj 

1 

I 

I 

I 

I 

( 

I 

I 



I 







e = .225, p = .0425 in, T = 70 F 




36 



V. 



CONCLUSION 



The moderate correlation among the constant load, con- 
stant displacement, constant displacement rate, and cyclic 
test cases tended to support the fracture postulate and 
gave credence to Hertzler's argument concerning dewetting* 
Test results, however, lacked absolute agreement as there 
was a 15% spread among the fracture characterizations for 
the four test cases. This disparity was primarily attri- 
buted to the nonlinearities inherent in the experimental 
tests. The effect of a finite crack tip radius on stress 
intensity factor was not known, but its influence did not 
ax>pear to be large since only a small range of crack tip 
radii were present in the experimental tests. 

The fracture postulate v/ould appear to be even more 
applicable to propellant, assuming the dewetting phenomenon 
could be incorporated into the fracture theory, because 
the relatively stiff nature of propellant and the con- 
straints imposed by the rocket motor structure would limit 
material strains within the linear assumptions used in 
the fracture postulate. 
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